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The solubility isotherm for CaHPCh ' 2H 2 (DCPD) in the three-component system Ca(OH) 2 
— H3PO4 — H 2 was determined in the pH range 3.5 to 6.8 by leaching a thermostated column of DCPD 
with dilute phosphoric acid solutions. In confirmatory experiments, equilibrium was approached both 
from super- and under-saturation by shaking DCPD with appropriate solutions. The calculated ionic 
activity product (Ca ++ ) X (HPQf), appeared to be a parabolic function of /;H with a minimum near pU 
5.0. The pH dependence of the ionic product could be accounted for by considering the ion pairs 
[CaHP04]° and [CaHiPCM* as semi-empirical parameters. Under the condition of saturation with 
respect to DCPD, the activity of the pair [CaHPQ 4 ]° must be a constant. The activity of the species 
[CaH 2 P04] + was shown to vary directly with hydrogen ion activity. The activities of the two ion pairs 
were adjusted to give a set ofpH-independent ionic activity products with a mean of 2. 19 ±0.01 1 X 10~ 7 . 
The stability constants for [CaHP0 4 ]° and [CaH 2 P0 4 ] + are 5.88 ±0.031 X 10 2 and 7.49 ±0.039, respec- 
tively. Experiments were conducted to study the hydrolysis of DCPD to more basic calcium phos- 
phates and the kinetics of these transformations is discussed. The significance of the ion pairs in 
human serum is considered. 

Key Words: Solubility of calcium monohydrogen phosphate dihydrate at 37.5 °C, calcium- 
phosphate ion-pairs, solubility of dicalcium phosphate dihydrate, solubility of dibasic 
calcium phosphate dihydrate. 



1. Introduction 

Dicalcium phosphate dihydrate (DCPD), CaHP0 4 
• 2H 2 0, has been found as a component of dental cal- 
culi [l] 1 and other biological deposits [2]. A role for 
DCPD as a nucleating agent [3] in calcification of 
tissues and its presence at the dental caries site [4] 
have also been suggested. Also, DCPD has chemical 
properties which are very helpful for the study of chem- 
ical interrelations among the various phosphates of 
biological significance; its high rate of solution for 
example, makes it possible to prepare reproducible 
calcium phosphate solutions within a wide range of 
concentrations and pH values which can be used in 
connection with crystal growth studies, calcification, 
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etc. It has commercial importance as a constituent 
of fertilizers, feed-grade dicalcium phosphate, and 
dentifrices. 

For these reasons, its solubility is of importance; 
there is disagreement, however, in reported values of 
its solubility product constant [5, 6, 7, 8, 9]. Perhaps 
the most important cause for inconsistent reports is 
its metastability with respect to the anhydrous salt, 
CaHPCX, and the more basic calcium phosphates, 
hydroxyapatite (HA), Cai (OH)2(PO 4 )6, and octacal- 
cium phosphate (OCP), Ca 8 H 2 (P0 4 )6 • 5H 2 0, [10]. 
Formation of complexes or ion pairs between the cal- 
cium and orthophosphate ions, as reported here, also 
contributes to uncertainty in solubility product con- 
stants calculated from solubility data. Its solubility 
isotherm in the three-component system, Ca(OH) 2 
— H3PO4— H 2 0, at 25 °C is known [9]. The solubility 
isotherm and solubility product at physiological tem- 
perature, 37.5 °C, are reported here. Equilibrations 
of DCPD with dilute phosphoric acids, using the leach- 
ing technique previously described [9], proved suc- 
cessful in minimizing formation of secondary phases. 
Batchwise equilibrations were used to verify the results 
obtained with the leaching technique. 
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2. Methods and Experimental Materials 
2.1. Preparation of Materials 

Dicalcium phosphate dihydrate was prepared by 
ammoniating in an ice bath an aqueous solution initially 
saturated with Ca(H 2 PQ 4 )2 ■ H 2 and CaHPQ 4 ■ 2H 2 0; 
details are described elsewhere [9]. Chemical anal- 
ysis of this salt yielded Ca, 23.20 percent and P, 17.93 
percent (theoretical: Ca, 23.20 percent and P, 18.02 
percent). Petrographically, the preparation appeared 
as pure, homogeneous DCPD. 

Phosphoric acid solutions were prepared from doubly 
crystallized hemihydrate, 2H3PO4 * H 2 0, and freshly 
boiled conductance water. The molarities of the final 
solutions were checked by chemical analysis. 

2.2. Leaching Experiments 

The leaching apparatus is shown in figure 1. An 
8-cm layer of DCPD crystals was poured on top of the 
fritted glass (medium porosity), F; phosphoric acid solu- 
tion was poured into the column; the solid was stirred 
to remove air bubbles; finally about 300 ml of phos- 
phoric acid solution was poured into the reservoir, S, 
continuosly flushing the apparatus with nitrogen. 
Water, from a constant temperature bath, was pumped 
through the jackets, W, maintaining the temperature 
at 37.5 ±0.1 °C as measured by thermometer T. The 
effluent from the column was collected in the thermo- 
stated cup D in which were mounted glass and calomel 
electrodes of a commercially available pH meter 
(claimed relative accuracy ±0.0037 pH units) and 
through which flowed a slow current of C0 2 -free and 
water-saturated nitrogen. When the pH of the solution 
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FIGURE 1. Leaching apparatus. 



entering the cup D became constant, several samples 
for analysis were collected through stopcocks C or E. 
The pH was rechecked several times during sampling 
and at various flow rates to be sure that the effluent was 
saturated. Standardization of the pH meter was made 
with certified NBS buffer standards [11]. The un- 
certainty in pH measurements was estimated to be 
about ±0.008 pH units. The hydraulic head was kept 
constant by means of the Mariotte valve M. The rate 
of flow at which the samples were collected was about 
25 ml per hour. When faster rates were needed, nitro- 
gen pressures of up to 3 psi were used to supplement 
the hydraulic head. 

2.3. End-Over-End Equilibrations 

A series of equilibrations was made in glass-stop- 
pered bottles, each of which contained 5 g of DCPD 
and 120 ml of phosphoric acid (1.00 X 10" 4 , 1.00 X 10" 3 , 
5.00 X10- 3 , or 5.00X 10" 2 M) by slowly rotating end- 
over-end in a water bath at 37.5 ±0.1 °C. In addition 
to the DCPD, 0.1 g of octacalcium phosphate, 
Ca 8 H 2 (P04)e ' 5H 2 0, (OCP) T was added to replicates of 
the two most dilute systems. At selected time inter- 
vals the contents of the bottles were filtered through the 
sintered glass in the thermostated apparatus shown in 
figure 1; pH values were determined and samples were 
collected for calcium and phosphorus analyses in the 
manner described above. 

Equilibrium was also approached from supersatura- 
tion, taking advantage of the negative thermal coef- 
ficient of solubility of DCPD: a column of DCPD in the 
apparatus shown in figure 1 was leached with 5.00 
X 10" 3 M phosphoric acid at 25 ±0.1 °C; thepH of the 
effluent was determined and a sample was taken for 
calcium and phosphorus analyses. The effluent was 
collected under a C0 2 -free atmosphere in a glass bottle 
containing about 2 g of DCPD; when the volume was 
about 120 ml, the bottle was stoppered and equilibrated 
in a water bath at 37.5 ±0.1 °C as described above. 
The pH and the calcium and phosphorus concentra- 
tions of the liquid phase were determined after 16 
days of equilibration. 

2.4. Calcium and Phosphorus Determinations 

Calcium was determined colorimetrically by the in- 
direct technique described by Banerjee et al. [12]. 
Phosphorus was determined using the vanadomolyb- 
date reagent of Brabson et al. [13]. In both analyses, 
samples were compared directly with blanks and stand- 
ards in a double beam spectrophotometer. Uncertain- 
ties in phosphorus and calcium determinations were 
estimated as 1.0 and 2.5 percent of the amounts 
present, respectively. 



3. Calculations 
3.1. Apparent Solubility Product Constants, K 8p 

The ionic strength, p,, is defined in these calculations 
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by 



lx = P[v + H + l2P+<t>/\] 



(1) 



in which P and H + represent molar concentrations of 
phosphorus and hydrogen ion and the other symbols 
are defined as: 



v =2MjP 



<S> 



(H+) 

'2K- 2 f f" 



4-4 



x=^+ ,H 



KiK> 



1+1 

K 2 f f" 



(2) 



(3) 



(4) 



In eqs (2), (3), and (4), (H + ) and M represent the 
hydrogen ion activity and the molar concentration of 
calcium, respectively, K\ and Kz are the first and sec- 
ond dissociation constants for phosphoric acid, and 
/' and /" are the molar activity coefficients for the 
species H 2 PQ 4 and HPQ7, respectively. The nu- 
merical values ior K\ and /Co used here were 6.53 X 10~ 3 
[14] and 6.57 XlO -8 , respectively. The value of Ki 
at 37.5 °C was interpolated by means of the equation 

pK 2 = 7. 1994 - 1.135X10- 2 X + 1. 998 XlO" 3 * 2 

-6.64X10- 5 X 3 with * = (T°K- 298. 15)/5 

derived from a portion (0-45 °C) of the experimental 
data of Bates and Acree [15]. 

Equations (1), (2), (3), and (4) are consistent with the 
assumptions that the only ionic phosphate species 
present in significant concentrations are H 2 PQ 4 and 
HPO| and that the contribution of hydroxyl ion to the 
ionic strength is negligible in comparison with that of 
the other ionic species present. These assumptions 
are justifiable in the pH range of the present investiga- 
tion. 

An approximate ionic strength was obtained by the 
use of eqs (1) to (4) assuming unitary values for/ ' and 
/", improved ionic activity coefficients, fi, were then 
calculated from the extension of the Debye-Hiickel 
limiting law, 



-\ogfr- 



AZfVfi 

1 + Ba,y/i* 



(5) 



The activity of HPO7 , (HPO4), was calculated from 

(HPO«) = />/X (6) 

and a solubility product for dicalcium phosphate di- 
hydrate, K f sp , was calculated according to eq (7), 
assuming a value of unity for the activity of water, 



K f 8p =Mf Ca (HPOj)=Mp4 £ - 

A 



(7) 



The solubility product calculated in this way was found 
to be dependent on concentration and is referred to 
here as the "apparent" solubility product. A more 
nearly constant value, K sp , was obtained by allowing 
for the presence of the ion pairs, [CaHP0 4 ]° and 
[CaH 2 P0 4 ] + as described below. 

3.2. Stability Constants 

The several stability and ionization constants con- 
sidered in this paper are defined as the equilibrium 
constants of the following reactions in solution: 

(a) Associative reactions (stability constants) 
Ca++ + HPO= <± CaHP0 4 [K*] (8) 
Ca++ + H 2 PO; *± Ca(H 2 P0 4 )+ [K*] . (9) 

(b) Dissociative reactions (ionization constants) 
Ca(H 2 P0 4 )+ *± CaHP0 4 + H+ [K y ] (10) 

CaHP0 4 *± CaPOj + H+ [K x ] . (11) 

For systems saturated with respect to DCPD, the 
concentration of the uncharged species [CaHP0 4 ]°, 
X°, is given by 

X° = K* (Ca + +)(HPO;)^. = K* KJU (12) 

in which (Ca ++ ) is the activity of ionic calcium and 
K sp is the activity product (Ca+ + )(HPOj). Therefore, 
X° is constant, assuming, as customary, that the ac- 
tivity coefficient, f x °, of an uncharged species in a 
dilute electrolyte solution is unity. The activity, 
(F + ), of the charged species [CaH 2 P0 4 ] + is given by 



in which Z\ represents the valence of the ith ion, a% its 
distance of closest approach, and A. and B are tem- 
perature-dependent constants; the numerical values 
chosen [16] for B and A were 0.331 X 10 8 and 0.522, 
respectively; the numerical values for a, were Ca ++ , 
6X10- 8 cm; HPOj and H 2 PO-, 4X10" 8 cm [17]. 
With activity coefficients obtained by the use of eq 
(5), a new value for /jl was calculated from eqs (1) to 
(4) and iteration was continued until convergence was 
attained, which usually occurred after three cycles. 



(Y+) = Y+ f v+ =X (M+)IKy 



(13) 



in which Y + represents the concentration of the species 
[CaH 2 P0 4 ] + ; accordingly, the activity of the charged 
species is proportional to the activity of hydrogen ion. 
In the calculations reported here, it was assumed that 
fy=f> It follows then, that if any two solutions 
designated by "a" and "i" are considered, the con- 
centration Yt in solution "6" is related to that in 
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solution "a" by the equation 

(H + ) & j^ 

b (m«f b - (14) 

For the calculations in this paper, solution "a" was 
taken as the most dilute solution in table 1. The 
selection of values of X° and Y% was made from an 
array of reasonable pairs set up for exploratory pur- 
poses; the choice of the values finally selected was 
based upon considerations outlined in the discussion. 



Table 1, Solubility of CaHP0 4 ■ 2H 2 0. Solution compositions 
obtained in leaching experiments 





Concentra- 


Composition of effluent 






System 

No. 






K' sp xW 


K sp xW 


H :i P0 4 










M x 10 4 


pU 


Ca, Mx 10 4 


P.Mx 10 4 






1 


5 


6.74k 


10.0, 


14.4« 


2.68 


2.15„ 


2 


10 


6.29.5 


13.7 4 


24.3,, 


2.59 fi 


2.22, 


3 


10 


6.29« 


13.9 7 


23.7 H 


2.59 4 


2.22. 


4 


50 


5.11, 


53.1 


101.4 


2.43« 


2.22 5 


5 


100 


4.58 7 


105.o 


203.-, 


2.41, 


2.14h 


6 


200 


4.12 5 


201.„ 


407.5 


2.57 2 


2.18 7 


7 


200 


4.13 5 


196.8 


403.5 


2.55« 


2.17 H 


8 


500 


3.56, 


484.2 


979.0 
S. D. a 

Avg. b 


2.865 
0.14 
2.59 ± 0.05o 


2.21, 
0.03 
2.19±0.01, 



a Standard deviation. 

b Averages given with the standard error of the mean. 



3.3. Corrected Solubility Product Constant, K sp 

The ionic strength was now calculated considering 
the balance equations (15) to (19) 



r = Y + +X° 
M = [Ca+ + ]-hr 

p = ^p(j)+ T 

[Ca ++ ] = M(l-/3) 
1 






pv 



(15) 
(16) 

(17) 

(18) 
(19) 



In eqs (15) to (19), r is the total concentration of cal- 
cium (or phosphorus) in ion pairs, [Ca ++ ] is the concen- 
tration of unpaired calcium; P® is the concentration 
of unpaired phosphate species of charge j; /3, the frac- 
tion of total calcium in paired forms, is f3 = r/M. The 
ionic strength, //,, is given by 



and the solubility product constant, K sp , is given by 

K sp = MP{\ -p) (1 - tfr) (f Ca lk). (21) 

For each selected pair of values for X° and Fj, eq 
(20) was solved by successive approximations starting 
with unitary activity coefficients analogously to the 
procedure outlined above. 

4. Results 

The composition of effluents in the leaching experi- 
ments and the corresponding calculated solubility 
products are shown in table 1. All the compositions 
are averages of at least two samples. Within experi- 
mental error, the final phosphorus concentration 
equaled the original phosphoric acid concentration 
plus the calcium concentration in the effluent, indi- 
cating that congruent dissolution took place in each 
case. 

The calcium and phosphorus concentrations in the 
various effluents at 37.5 °C are plotted against pH in 
figure 2; corresponding curves for 25 °C [9} are shown 
as broken lines. It is apparent that DCPD is more 
soluble at 25 °C than at 37.5 °C. 

Apparent solubility products, Kg P , listed in table 1, 
are plotted against pH in figure 3. The concave shape 
of this curve is consistent with a model in which the 
iort pairs [CaHP0 4 ]° and [CaH 2 P04] + are present in 
solution. The corrected solubility product constants, 
K sp , calculated on the basis of this model are also 
plotted in figure 3. The maximum variation in K sp , 
about 4 percent, is satisfactory in view of the analytical 
uncertainties. In contrast, Kg P varied by about 20 
percent. 

In table 2 are listed the parameters used in the cal- 
culations along with pertinent stability and ionization 



M = P 


_»-§g*) + v»A)(l-!*) 




-(1/2P)(J°-H + ) 


(20) 
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FIGURE 2. Solubility isotherm of DCPD at 37.5 °C,full 
line, and at 25 °C, broken line [9]. 

Open circles, compositions obtained in leaching apparatus; triangles, 
compositions obtained in end-over-end experiments; squares, initial and final 
compositions in super-saturation experiments (see tabulation in text). 
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FIGURE 3. Solubility product constants (apparent and 
corrected) and percent P in ion pairs plotted against 
pH. 

Calculated errors, based on analytical uncertainties, are shown for each 
point. 



Table 2. Ionic strengths and concentrations and stability constants 
of ion pairs a 



System 


fix 10* 


tX 10 4 


(8 


Y + x 10 4 


K* x 10- 2 


K* 


No. 














1 


0.300 


1.34, 


0.134 


0.0475 


6.00,, 


7.64, 


2 


.408 


1.42,; 


.104 


.133, 


5.80:, 


7.39, 


3 


.410 


1.42, 


.102 


.136, 


5.80 (i 


7.39 4 


4 


1.51,, 


3.49, 


.066 


2.20 8 


5.79« 


7.38., 


5 


2.94, 


8.95, 


.085 


7.65h 


6.00« 


7.64» 


6 


5.57, 


24.5* 


.120 


23.2* 


5.89 4 


7.51a 


7 


5.46 4 


24.0, 


.124 


22.7, 


5.92 H 


7.55, 


8 


12.6 () 


92.9« 


.192 


91.6h 


5.83., 


7.43, 



S. D. 1 
Avg. c 



0.09 0. 1 1 

5.88 ±0.03, 7.49 + 0.03* 



a Values in table calculated for the final set X°= 1.293 X 10~ 4 and Y+ = 0.0475 X 10" 4 M 

''Standard deviation. 

'Averages given with the standard error of the mean. 

The value for K u was 5.160 X lQ-«. 





Table 3. 


End-over- 


?nd equi 


librations 










Solution composition 




Concentra- 
tion of 


Equilibration 

time — hours 






phases 










II.I'O, 




pH 


Ca. Mx 10 4 


/>, MX 10 4 




Mx 10 2 


264 


5. OK, 


19.8, 




DCPD+OCP 


0.01 


34.9, 


DCPD + OCP 


.1 


288 


5.65j 


28.1a 


52.7, 


DCPD 


.01 


5 


7.42 H 


7.20 


8.24 


DCPD 


.01 


24 


^.HH, 


20.2k 


37.9< 


DCPD 


.01 


48 


5.80,, 


22.7 a 


42.6, 


DC1M) 


.01 


240 


5.55a 


30.5 7 


58.4 4 


DCPD 


.01 


000 


5.28, 


43.2, 


84.8, 


DCPD 


0.5 


192 


5.13, 


51.4„ 


99.9 4 


DCPD 


5.0 


192 


3.57„ 


473.H 


963.„ 



two of the systems, the last two entries in table 3, 
showed congruent dissolution. These contained only 
DCPD and were equilibrated with 5.00X10- 3 M and 
5.00XlO _2 /lf phosphoric acid, respectively. The 
compositions of these two solutions are the only ones 
that are essentially the same as those obtained in the 
corresponding leaching experiments, table 1. 

In the OCP-free systems equilibrated with 1.00 
X 10" 4 M phosphoric acid, table 3, a more basic phos- 
phate precipitated as evidenced by the changes in com- 
position with time. The rate of this precipitation was 
relatively fast at the beginning and then decreased. 

The two systems containing both DCPD and OCP, 
table 3, were included with the hope of obtaining the 
singular point for the two salts (i.e., the solution which 
would be in equilibrium with both salts). The lack 
of agreement in the final composition of the two solu- 
tions indicates that at least one of them did not remain 
at the singular point composition. Hydrolysis of OCP 
to a more basic phosphate, presumably hydroxyapatite, 
would account for this result. 

The initial and final compositions in the experiment 
in which equilibrium was approached from supersatu- 
ration are shown in the tabulation: 



constants. The values for K% and K y were calcu- 
lated using eqs (12) and (13), respectively. The values 



g* _ (Y + )K 2 



(22) 



The constancy in the value of K y is imposed by the 
calculative procedure, as will be explained later. 
The near constancy in the values of K* and K^ (table 2) 
lends strong support for the existence of the complexes 
postulated here. 

The fraction of phosphorus complexed, t/P, is also 
plotted in figure 3. There appears to be a minimum 
of complexing at pH 5. As might be anticipated, the 
minimum in the fraction of P present as ion pairs is 
nearly coincident with the minimum in the K' 8p curve 
in the same figure. 

The compositions of all the solutions obtained in the 
end-over-end equilibrations, listed in table 3 and shown 
in figure 2 as triangles, lie close to the isotherm ob- 
tained in the leaching experiments. However, only 



Initial 
Final 



pH 



5.20 2 
5.14a 



Ca 
MX 10 4 



52.9 8 
48.2 2 



P 

Mx 1() 4 



101. 9 
100.8 



The initial and final compositions are shown as squares 
in figure 2. It is apparent that the former lies on the 
25 °C isotherm and the latter on the 37.5 °C isotherm. 



5. Discussion 

5.1. Leaching Experiments 

The more basic calcium phosphates (such as HA and 
OCP) were not formed in the leaching experiments as 
long as the phosphoric acid solutions were more con- 
centrated than 1.00X10~ 4 M. Attempts to attain 
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equilibrium with 1.00X10~ 4 M acid failed because 
considerable hydrolysis of DCPD was evident even at 
flow rates 3 to 4 times higher than the rate used with 
more concentrated acids. This behavior is in marked 
contrast to that at 25 °C [9] which showed no evidence 
of hydrolysis with the same acid. 

The following facts show that equilibrium was 
closely approached in the experiments reported here: 
(a) the plots of experimental points yielded smooth 
curves; (b) the end-over-end experiments (systems 8 
and 9, table 3), which represent much longer equilibra- 
tion times, gave solution compositions in agreement 
with those obtained in the leaching experiments; 
(c) the system in which equilibrium was approached 
from supersaturation, agreed with that in which equi- 
librium was attained from undersaturation with the 
same acid; (d) within experimental error, a constant 
solubility product was obtained when corrections for 
the presence of complexes in solution were introduced. 

Dicalcium phosphate dihydrate is metastable with 
respect to anhydrous dicalcium phosphate over the 
entire range studied and with respect to OCP and HA, 
in part of the range. The material in the column was 
examined after each experiment with the petrographic 
microscope. With one exception, the system in which 
the most concentrated acid was used, no anhydrous 
dicalcium phosphate was found. The trace found 
in the one exception is believed to be too small to 
influence the corresponding point on the isotherm 
because it is known that the rate of precipitation of 
the anhydrous salt is much less than the rate of dis- 
solution of DCPD [18]. 



5.2. End-Over-End Experiments 

Since all the compositions shown in table 3 fall on 
the DCPD isotherm (fig. 2) obtained in the leaching 
experiments, the rate of precipitation of more basic 
calcium phosphates (such as OCP and HA) in the 
end-over-end experiments must have been signifi- 
cantly lower than the rate of dissolution of DCPD. 

With only DCPD and OCP present as solid phases, 
the solution composition should approach the DCPD- 
OCP singular point. The formation of HA in such a 
system would cause compositions to shift toward the 
more acid region. This would create a condition of 
undersaturation with respect to OCP; the OCP would 
then dissolve, leading to supersaturation with respect 
to DCPD. The dissolution of OCP and precipitation 
of DCPD would thus tend to maintain the pH at a 
higher level, approaching that of the OCP-DCPD 
singular point. Therefore, depending upon the rela- 
tive rates of these three reactions, the solution com- 
position would be either at the DCPD-OCP singular 
point or somewhere on the acid side of this point. In 
agreement with these considerations, the 1.00 X 10~ 3 M 
acid yielded an equilibrium pH of 6.295 in the leaching 
experiments; the same acid after 288 hr in the end-over- 
end experiments gave a pH of 5.653. The difference 
must be attributed to the formation of HA. Likewise, 
it is probable that the DCPD-OCP singular point has 



a pH above the value 5.984 obtained with the 1.00 
X10- 4 M acid after 264 hr. 

The ionic product for OCP, (Ca+ + ) 4 (H+)(POf) 3 , 
calculated for the solution with pH 5.984 (/3=0.08) 
is 3,8 X 10 -48 . Because of the foregoing reasons, the 
solubility product of OCP at 37.5 °C must be greater 
than 3.8 X 10~ 48 which puts it in the neighborhood of 
the product reported for 25 °C, 1.25 X lO 47 [10]. 

In the absence of OCP, the 1.00 X 10~ 4 M acid gave 
solutions with pH values below 5.984, the value 
obtained in the presence of OCP; at 24 hr thepH was 
5.883 and it decreased to 5.281 after 960 hr. Pre- 
cipitation of HA liberates acid which is not neutralized 
by dissolution of OCP when the latter is absent. 
Consequently, the pH in the time series with 1.00 
X 10 -4 M acid in the absence of OCP, decreased fairly 
rapidly to approach that of the DCPD-HA singular 
point. The results suggest that the pH of the DCPD- 
HA singular point is below 5.281, and probably above 
5.132, the value obtained in system 8, table 3, in which 
no hydrolysis was detected after 192 hr. 

The hydrolysis steps outlined in the foregoing dis- 
cussion are consistent with the known behavior of 
DCPD [9, 10] and of OCP [19]. The formation of 
OCP by hydrolysis of DCPD in aqueous solutions 
occurs at 25 °C, but it was possible [10] to determine 
the singular point for these two salts, because at 25 °C 
hydrolysis of OCP to HA is too slow to displace the 
composition from that of the singular point. 



5.3. Ion-Pair Formation 

The values for the concentrations of the ion-pairs 
[CaHP0 4 ]° and [CaH 2 P0 4 ] + , *° and F + , in the most 
dilute solution were selected from an array according 
to the following criteria: 

(a) The calculated value for K sp cannot be higher 
than the minimum value for K' 8p in table 1; (b) the sum 
of the concentrations of the ion pairs cannot be higher 
than either the total Ca or P concentrations in any of 
the systems considered; (c) there must be a minimum 
in the fraction (r/P) of P complexed and it should coin- 
cide with the minimum in the value of KL in figure 3; 
(d) the variance of K gp obtained with the final values 
of X° and Y+ should be an absolute minimum with 
respect to all other possible pairs of values for these 
parameters. 

Designating values for X° and Y + for the most dilute 
solution fixes a numerical value of K y for all the other 
solutions because the quantities (Y + )/(H + ) and X° in 
eq (13) are constants in the adopted model. This 
is why the quantity K y does not appear in table 2 along 
with the calculated values for K% and Ky. The value 
derived for K y is 5.16 X 10~ 6 . 

The value for the acid ionization constant, K y , of 
the species [CaH 2 P0 4 ] + , 5.16 X 10~ 6 , is one hundred 
times the second ionization constant of phosphoric 
acid; this ion-pair is thus a stronger acid than H 2 PO^, 
presumably due to the proximity of the charge of the 
calcium ion. 
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If the reasonable assumption is made that the ratio 
of KyjKx is of the same order of magnitude as KJK>> 
or K2IK3 for phosphoric acid, the value of K x would be 
of the order of 10 -11 . This means that the activity 
of the species [CaP0 4 ]~ in the system with the highest 
pH would be of the order of 10~ 8 M, which is negli- 
gible relative to the activities of Ca + + and HPOf. 
Therefore, ignoring the concentration of the species 
[CaP0 4 ]~ in the calculations is justified. 

In a previous study of the solubility of DCPD at 
25 °C [9], the possibility of complexing was not ruled 
out, but the data were treated by an empirical method 
that gave K sp by extrapolation to zero ionic strength. 
The 25° data are being re-evaluated in this laboratory, 
taking into account ion-pair formation. 

According to the present treatment, it appears that 
as much as 20 percent of the calcium was ion-paired 
(e.g., system No. 8, table 1). The model adopted for 
the treatment of the ion pairs is not exempt from un- 
certainties; for example, it is impossible to ascertain 
the adequacy of eq (5) for calculating activity coeffi- 
cients. On the other hand, the variation of the solu- 
bility product constant could not be explained simply 
in terms of activity coefficients because (a) the ex- 
tended Debye-Hiickel limiting law should give reason- 
ably good values for activity coefficients in solutions 
as dilute as those used here, and (b) it is unlikely that 
the minimum in the value of K' sp which occurs at 
jjl ~ 0.025 is related to minima in mean activity coeffi- 
cients which are sometimes observed; such minima 
occur at concentrations much higher than those 
described here. 



5.4. Comparison With Reported Data 

Davies and Hoyle [20] reported dissociation constants 
at 25 °C for systems containing sodium and potassium 
phosphate buffers and calcium iodate. Allowing 
for experimental differences, the constants obtained 
by Davies and Hoyle, £* (25°)= 1/.002 = 5X 10 2 and 
Ky (25°)= 1/. 084 =11.8, are in reasonable agreement 
with those given here, K* x (37.5°) = 5.88 X 10 2 and K* 
(37.5°) = 7.49. 

Bjerrum [21] studied the three component system 
Ca(OH) 2 — H3PO4— H,0 at 18 °C in the pH range 6 to 8 
in the absence of solid phases. He concluded that a 
model involving the ion-pair [CaHP0 4 ]° or any of the 
ion pairs [Ca(HP0 4 ) 2 ] 2 -, [Ca(HP0 4 )(P0 4 )] 3 -, or 
[Ca(P0 4 ) 2 ] 4 ~, would fit his data. Models involving 
[CaHP0 4 ]° alone or in combination with any of the 
above ion pairs with Ca/P ratios of 1/2 gave variable 
K sp values when used with the present data. The 
only model giving reasonably constant K sp values was 
the one with the ion pairs [CaHP0 4 ]° and [CaH 2 P0 4 ] + . 
It is believed that the condition of saturation with 
respect to DCPD, which demands a constant activity 
for the species [CaHP0 4 ]°, is more determinative 
than the conditions used by Bjerrum. 

More recently, Greenwald [22J used Bjerrum's data 
to calculate values for the instability constant of 



[CaHP0 4 ]°; these varied from 1.01 to 9.55 X 10" 3 . As 
in previous reports by the same author [23, 24], ionic 
concentrations were used instead of activities. For 
this reason, and because the procedure used in his 
calculations is not detailed, a direct comparison of 
his results with those given here is not possible. His 
average value for the instability constant, 2.77 X 10 -3 , 
gives a K* (18 °C)= 1/(2.77 X 10 : *) = 3.61 X 10 2 . This, 
too, gives order of magnitude agreement with the 
value of K* (37.5 °C) = 5.88 X 10 2 given here. 

Similar rough agreement is found with the result 
of Gosselin and Coghlan [25], K* (25 °C) = 1/1.38 X 10" 2 
= 0.72 X 10 2 , who also used ionic concentrations in- 
stead of activities in their calculations. 

The negative temperature coefficient for the solu- 
bility of DCPD relates to properties of both the solu- 
tion and the solid. For example, with an increase in 
temperature the activity of water increases and it is 
conceivable that less water is associated with the 
calcium and phosphate ions which, then, for com- 
parable concentrations, would display a higher activ- 
ity; this should result in lower solubility of DCPD as 
the temperature increases. However, if the solubility 
product, Kgp 9 decreases as the temperature increases, 
the magnitude of the standard free energy of formation 
of DCPD must increase with temperature more rapidly 
than the sum of the free energies of formation of 
Ca+ + , HPOf, and 2H 2 0. The solubility product of 
DCPD at 25 °C [9], 2.77 X 10" 7 , is higher than the value, 
2.19X 10~ 7 , at 37.5 °C. However, since no consider- 
ation was given to the formation of ion pairs in solution 
at 25 °C, a meaningful comparison is not possible at 
present. 

An important application of the stability constants 
relates to the formation of ion pairs in biological 
fluids. Representative figures cited for ionic Ca ++ 
and HPOj activities in adult human serum, 4.7 X 10 ~ 4 
and 1.9 X 10~ 4 M, respectively [26), give the result at 
pH = 7.4 that about five percent of the calcium or 
inorganic phosphorus is in the form of the ion pair 
[CaHPQ 4 ]°. In infant serum, with Ca + + and HPQf 
activities of 4.7 X 10~ 4 and 3.7 X 10~ 4 , the percentage 
of calcium and phosphorus in paired forms would 
amount to 8 and 5, respectively. For both adults 
and infants, the concentrations of the ion pairs 
[CaH 2 P0 4 ] + and [CaP0 4 ]" in serum appear to be 
negligible. 
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